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A search for production of supersymmetric particles in final states containing jets, missing trans- 
verse momentum, and at least one hadronically decaying r lepton is presented. The data were 
recorded by the ATLAS experiment in y/s = 7 TeV proton-proton collisions at the Large Hadron Col- 
lider. No excess above the Standard Model background expectation was observed in 2.05 fb _1 of 
data. The results are interpreted in the context of gauge mediated supersymmetry breaking models 
with M mess = 250 TeV, N 5 = 3, fi > 0, and C grav = 1. The production of supersymmetric particles is 
excluded at 95% C.L. up to a supersymmetry breaking scale A = 30 TeV, independent of tan/3, and 
up to A = 43 TeV for large tan/3. 
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1. Introduction 

Supersymmetry (SUSY) \Mk 

is a well-motivated theo- 
retical concept that introduces a symmetry between bosons 
and fermions. As a consequence, every Standard Model 
(SM) particle has a SUSY partner with the same mass 
and quantum numbers except for the spin which differs 
by half a unit. Since none of these partners has been ob- 
served SUSY must be a broken symmetry if realised in 
nature. If R-parity is conserved I3~14|, SUSY particles 
can only be produced in pairs and would decay through 
cascades involving lighter SUSY particles. These decay 
cascades end in the production of the lightest supersym- 
metric particle (LSP), which is stable and escapes the de- 
tector unseen, giving rise to missing transverse momentum 
in the detector. SUSY can remedy various shortcomings 
of the Standard Model, such as the hierarchy problem [lj- 
[lij ]. the lack of a dark matter candidate 20. 21 1 and the 
non- unification of the gauge couplings 22h25| . To achieve 
this, the masses of at least some SUSY particles must be 
near the weak scale, and therefore, if weak-scale SUSY is 
realised in nature, there are good prospects to discover it 
at the Large Hadron Collider (LHC). 

In certain SUSY models, large mixing between left and 
right sfermions, the partners of the left-handed and right- 
handed SM fermions, implies that the lightest sfermions 
belong to the third generation. This leads to a large pro- 
duction rate of r leptons from decays of f sleptons and 
gauginos, the partners of the SM gauge bosons, in SUSY 
cascade decays. For example, in the context of Gauge Me- 
diated SUSY Breaking (GMSB) [H-|3l| the lighter of the 
two f sleptons is the next-to-lightest supersymmetric par- 
ticle (NLSP) for a large part of the parameter space, and 



the very light gravitino, G, is the LSP. Hence f sleptons 
decay to a t lepton and a gravitino. While this f — > tG 
process is the dominant source of r leptons from SUSY de- 
cays in certain regions of GMSB model parameter space, 
the analysis presented here is sensitive to any process pro- 
ducing r leptons in association with jets and missing trans- 
verse momentum. 

This Letter presents a search for supersymmetry in fi- 
nal states with at least one hadronically decaying r lepton, 
missing transverse momentum and jets with the ATLAS 
detector at the LHC. The results of the search are inter- 
preted within the GMSB model. Previous experiments at 
LEP [32 - 34 1 have placed constraints on f and e masses and 
on more generic GMSB signatures. Among these the limits 
from the OPAL experiment 32j were the most stringent, 
excluding f NLSPs with masses below 87.4 GeV. The DO 
Collaboration performed a search for squark production 
in events with hadronically decaying r leptons, jets, and 
missing transverse momentum [35j, and the CMS Collab- 
oration performed searches for new physics in same-sign 
ditau events [3(| and multi-lepton events [3lJ including r 
pairs, but the GMSB model was not specifically considered 
in any of these results. A search for supersymmetry in fi- 
nal states containing at least two hadronically decaying r 
leptons, missing transverse momentum, and jets with the 
ATLAS detector is presented in another Letter 



2. ATLAS detector 

The ATLAS detector [39j is a multipurpose particle 
physics apparatus with a forward-backward symmetric cylin- 
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drical geometry and nearly Att coverage in solid angle Q. 
The inner tracking detector consists of a silicon pixel de- 
tector, a silicon microstrip detector, and a transition ra- 
diation tracker. The inner detector is surrounded by a 
thin superconducting solenoid providing a 2 T axial mag- 
netic field and by high-granularity liquid-argon sampling 
calorimeters. An iron-scintillator tile calorimeter provides 
hadronic coverage in the central rapidity range. A muon 
spectrometer consisting of large superconducting toroids 
and a system of precision tracking chambers surrounds the 
calorimeters. 

3. Data and simulated samples 

The analysis is based on data collected by the ATLAS 
detector in proton-proton collisions at a centre-of-mass en- 
ergy of 7 TeV between March and August 2011. Applica- 
tion of beam, detector, and data-quality requirements re- 
sulted in an integrated luminosity of 2.05 ± 0.08 fb -1 [Ijl 
4l|. The data were collected using triggers based on one 
jet with transverse momentum p^ > 75 GeV, measured 
at the raw electromagnetic scale, and missing transverse 
momentum above 45 GeV. 

In GMSB models, the breaking of SUSY is mediated 
through flavour-blind SM gauge interactions of messenger 
fields with mass scale M mcss which is small compared to 
the Planck mass. In addition to M mess , the free parame- 
ters in GMSB models are the scale of the SUSY breaking, 
A, the number of messenger fields, N5, the sign of the 
Higgsino mixing parameter, sign(^), the scale factor for 
the gravitino mass, C grav , and the ratio of the vacuum 
expectation values of the two Higgs doublets, tan /3. In 
this analysis, GMSB models are studied in the A — tan /3 
plane for fixed M mcss = 250 TeV, N 5 = 3, sign(/x) = +1 
and Cg ra v = 1- The chosen set of parameter values re- 
stricts the analysis to specific final states relevant for the 
search with r leptons and to promptly decaying NLSPs. 
For N$ > 2 and large tan/? the lightest f slepton, fi, is 
the NLSP. 

Samples of simulated GMSB events are generated with 
the Herwig++ [42| generator for ten values of A in the 
range 10 < A < 85 TeV and ten values of tan /3 in the 
range 2 < tan/3 < 45, with the SUSY mass spectra gen- 



erated using ISA JET 7.80 [43(. The MRST2007 L0* [44 1 
parton distribution functions (PDFs) are used. The pro- 
duction cross sections are calculated with PRDSPIND [45- 
|48] | to next-to-leading order in the QCD coupling using 
the next-to-leading-order CTEQ6 . 6 49] PDF set. The two 



1 ATLAS uses a right-handed coordinate system with its origin 
at the nominal interaction point in the centre of the detector and 
the 2-axis coinciding with the axis of the beam pipe. The rr-axis 
points from the interaction point to the centre of the LHC ring, and 
the y-axis points upward. Cylindrical coordinates (r, rfi) are used in 
the transverse plane, <f> being the azimuthal angle around the beam 
pipe. The pseudorapidity is defined in terms of the polar angle 9 as 
r\ = -lntan(0/2) 



samples with A = 30 (40) TeV and tan/3 = 20 (30), which 
have cross sections of 1.95 (0.41) pb, are used as represen- 
tative points for the optimization of the event selection. 

The dominant background processes in this search are 
production of W and Z bosons in association with jets 
(W+jets and Z + jets), top quark pair (tt) and single top 
quark production. The W + jets and Z + jets production 
processes are simulated with the ALPGEN [H3| generator, 
using the CTEQ6L1 [5l| PDF set, and are normalized to a 
cross section of 31.4 nb and 9.02 nb 52-54j], respectively. 
The ti, single-top and diboson production processes are 
generated with MC0NL0 [55| and the CTEQ6.6 [H PDF 
set, and are normalized using a cross section of 0.165 nb, 
0.085 nb [H-[5i] and 0.071 nb respectively. Parton 

showers and hadronization are simulated with HERWIG and 
the underlying event is modelled with JIMMY [6l| . The pro- 
grams TAUDLA [H EH and PHOTOS [H are used to model 
the decays of r leptons and the radiation of photons, re- 
spectively. T he p roduction of multijet events is simulated 
with PYTHIA 65J, though the multijet background yield in 
this analysis is estimated using data. All simulated sam- 
ples are processed through a full simulation of the ATLAS 
detector 66] based on Geant4 67]. To match the pile-up 
(overlap of several interactions in the same bunch crossing) 
observed in the data, the generated signal and background 
events are overlaid with minimum-bias events [68], |69[ and 
the resulting events are reweighted so that the distribution 
of the number of interactions per bunch crossing agrees 
with the data. 

4. Object reconstruction 

Jet candidates are reconstructed with the anti-fct clus- 
tering algorithm (70I ] with radius parameter R — 0.4. The 
inputs to this algorithm are clusters of calorimeter cells 
seeded by cells with energy significantly above the mea- 
sured noise. Jets are constructed by performing a four- 
vector sum over these clusters, treating each cluster as a 
four- vector with zero mass. Jets are corrected for calorime- 
ter non-compensation, upstream material, and other ef- 
fects using pY- and ^-dependent correction factors obtained 
from Monte Carlo simulation and validated with exten- 
sive test-beam and collision-data studies 7l|. Only jet 
candidates with pt > 30 GeV, \rj\ < 2.8 and a distance 
AR > 0.2 with respect to the nearest identified electron 
are considered as real hadronic jets, where the distance is 
defined as AR = ^(Arj) 2 + (Ac/)) 2 . 

The electron and muon identification criteria are iden- 
tical to those in Ref. [72j. Electrons and muons are only 
considered if they satisfy pt > 20 GeV and AR > 0.4 with 
respect to the nearest identified jet. 

The magnitude of the missing transverse momentum, 
E^ lss , is computed from the vector sum of the transverse 
momenta of all identified electrons and muons, all jets, and 
remaining clusters of calorimeter cells with \rj\ < 4.5 ;73|. 

Hadronically decaying r leptons are reconstructed from 
jet candidates with pt > 10 GeV and are distinguished 
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from quark- or gluon- initiated jets using a boosted decision 
tree (BDT) based on eleven discriminating shower-shape 
and tracking variables [74j ]. Electrons are further rejected 
using transition radiation and calorimetric information. 
An energy calibration factor for hadronically decaying r 
leptons is applied as function of px and rj. Candidates are 
required to satisfy p T > 20 GeV and \rj\ < 2.5 and to have 
one or three associated reconstructed tracks (prongs) with 
total charge ±1. The r candidates are required to sat- 
isfy a px-dependent BDT output criterion [74f chosen to 
give ~ 30% (~ 50%) signal efficiency for one-prong (three- 
prong) r candidates as estimated in Z(—¥ rr) +jets events. 
The BDT selection has a corresponding background ac- 
ceptance of ~ 0.5% (~ 3%), estimated in dijet events, and 
the different selection criteria reflect different abundances 
of one- and three-prong jets in background samples. 

During a part of the data-taking period, an electron- 
ics failure in the liquid-argon calorimeter created a dead 
region in the second and third layer of the calorimeter, 
corresponding to approximately 1.4 x 0.2 rad in Ary x A(f>. 
A correction is made to the jet energy using energy deposi- 
tions in cells neighbouring the dead region; events having 
at least one jet, including the leading r candidate, in this 
region for which the corrected energy is above 30 GeV are 
discarded, resulting in a loss of ~ 6% of the data sample. 

5. Event selection 

Events are required to have a reconstructed primary 
vertex with at least five associated tracks with px > 500 MeV. 
Events are rejected if they contain identified electrons or 
muons or if any jet or r candidate is consistent with aris- 
ing from detector noise or non-collision background (7lj . 
Events are required to contain one or more identified r 
candidates, at least two jets, one with px > 30 GeV and 
another with px > 130 GeV, and missing transverse mo- 
mentum E™ ss > 130 GeV. The latter two requirements 
ensure that the trigger efficiency is above 98% in both data 
and simulation. 

The two jets leading in px are required to be sepa- 
rated in azimuth from the direction of the missing trans- 
verse momentum by more than 0.3 rad. This requirement 
reduces multijet events, which typically have instrumen- 
tal missing transverse momentum aligned with the leading 
jets. Multijet events are further suppressed by requiring 
E^ lss /m c ff > 0.25, where the effective mass, m c g, is de- 
fined as the scalar sum of E™ lss , the px of the two leading 
jets, and the p T of the leading r candidate. 

Events are required to have a transverse mass, mx, 
above 110 GeV. The transverse mass is defined as 

m T = J ml + 2p T £" iss (l - cos A0(p T , E™ iss )) , 

where A0(p T , E™ lsa ) is the azimuthal angle between the 
r and the direction of the missing transverse momentum. 
This requirement suppresses backgrounds due to VK+jets 
and top-quark production. The remaining SM backgrounds 



are further suppressed by requiring m e s > 600 GeV. This 
is the final selection defining the signal region for the anal- 
ysis. The mx and m e s requirements as well as the criteria 
used for the suppression of multijet events are chosen to 
maximise the signal significance computed with the Asi- 
mov approximation |75| . 

6. Background estimation 

Background processes are divided into three classes 
which are estimated separately: events with true r leptons 
from t —¥ brv decays (both top-quark-pair and single top 
quark production) and W{-) tv) + jets events; events with 
misidentified ('fake') r candidates in top, W + jets, and 
Z + jets events; and events with fake r candidates in mul- 
tijet events. The two fake-r classes are treated separately 
to account for differences in r misidentification probabili- 
ties due to different event topologies and jet composition. 

Events with true r leptons are estimated in a con- 
trol region defined by replacing the requirement on the 
transverse mass in the final selection with the requirement 
mx < 70 GeV. For events with a correctly reconstructed 
r lepton and with E™ lss entirely due to a single neutrino, 
mx is kinematically bounded from above by the W mass, 
within the detector resolution; by requiring mx < 70 GeV, 
more than 90% of the events in the resulting control re- 
gion are expected to contain true r leptons from top-quark 
and W decays. The composition of the event sample in 
this control region is given in Table [TJ Within this con- 
trol region, the background due to Z decays is estimated 
from simulation and the remaining small background due 
to multijet events is estimated using a procedure similar 
to that used to estimate the multijet background in the 
signal region, described below. 

Within the mx < 70 GeV control region, top-quark 
and VF+jets yields are estimated individually with a maximum- 
likelihood fit to the output distribution of a BDT built 
from four variables: the number of 6-quark jets, the total 
jet multiplicity, the transverse momentum of the second- 
leading jet, and the transverse thrust T of the event, de- 
fined as T = maxftl^n -PT,i/\ 2jFr,i|}, where i runs 
over the missing transverse momentum and all jets, ex- 
cluding the tau candidates, with transverse momentum 
vectors px,i, and the transverse thrust axis is given by 
the unit vector n for which the maximum is attained. 
Top-quark events have more reconstructed 6-quark jets, 
a higher jet multiplicity, higher jet momenta, and tend to 
be more spherical than W + jets events. Jets containing b 
quarks arc identified with about 60% efficiency, evaluated 
with top-quark events, using secondary vertex reconstruc- 
tion and three-dimensional impact parameters of tracks 
associated with the jet j76[. The output distribution of 
this BDT is shown in Figured] along with the results of the 
fit. The results of the fit are scale factors for W-^+jets and 
top quark backgrounds which reflect differences in cross 
sections and reconstruction efficiencies between data and 
simulation. The measured scale factors are 1.22 ± 0.13 for 
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top events and 0.71 ±0.03 for TV + jets events. These scale 
factors are applied to simulated event samples in the sig- 
nal region to derive the final expected true-r yields from 
background processes. 

Table 1: Numbers of observed and expected events in the true-T- 
dominated W/top control region, defined as mx < 70 GeV. The 
numbers shown for W+jets and top are from Monte Carlo simulation 
and do not include the correction factors derived from this control 
region. The correction factors obtained from a fit to data are 1.22 ± 
0.13 for top and 0.71 ± 0.03 for W + jets. The true-r purity is 97% 
for top, 96% for W + jets and 87% for Z + jets. 
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Data 
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BDT output 

Figure 1: Output distribution of the BDT used to discriminate 
W + jets from top-quark events in the low-mx control region, de- 
fined as rri'Y < 70 GeV. Background distributions are taken from 
simulation. The yield for W + jets and top backgrounds are taken 
from a maximum-likelihood fit to this distribution. The solid (red) 
line with shaded (yellow) error band corresponds to the total SM 
prediction, while the points are data. 

For the estimation of backgrounds due to fake r can- 
didates in top-quark, W + jets, and Z + jets events, a 
second control sample is defined by selecting events that 
fulfill the event selection but with modified criteria on rm 
and m c g: itlt > 70 GeV and either mx < 110 GeV or 
m e ff < 600 GeV. Since the m T distribution falls off rapidly 
above the W mass for true-r events, the intermediate jtit 
region selected here is relatively enhanced in fake-r events, 
and the overall composition of this region is expected to be 
very similar to that of the signal region. Multijet events 
are expected to make up less than 3% of this sample and 
are estimated from simulation. The composition of the 
fake-r-enhanced sample in this control region is shown in 
Table [2] Within this control region, true-r backgrounds 



are subtracted using estimates derived from the true-r- 
dominated control region. The numbers of events remain- 
ing after the true-r subtraction are used to determine a 
scale factor, 0.50±0.08, which is then applied to simulated 
samples of fake-r events in the signal region to obtain a 
final background estimate. While this scale factor differs 
significantly from unity, it is consistent with other ATLAS 
studies of the performance of r fake rates in simulation. 

Table 2: Numbers of observed and expected events in the fake-r- 
enhanced control region. The numbers of expected W + jets and 
top-quark events have been corrected by the correction factors mea- 
sured in the true-r-dominated region. The fake-r correction factor 
obtained from data is 0.50 ± 0.08. 





True r 


Fake r 


Total 


Top 


53.3 ± 7.5 


37.8 ± 5.8 


91.1 ± 9.4 


W + jets 


80.5 ± 6.9 


33.3 ± 4.1 


113.8 ± 8.0 


Z + jets 


5.1 ± 1.6 


41.5 ± 10.8 


46.6 ± 10.9 


Multijet 


± 


2.9 ± 1.0 


2.9 ± 1.0 


Total 


139 ± 10 


116 ± 13 


254 ± 17 


Data 


197 



Backgrounds due to multijet events are estimated in a 
third control region in which either E^ lss /m c s < 0.25 or 
one of the two leading jets is aligned in azimuth with the 
missing transverse momentum direction. Within this sam- 
ple, the probability for jets (which contain very few true 
r leptons) to satisfy the r selection criteria is estimated 
by applying the selection to randomly chosen jet candi- 
dates. This probability is then applied to a complemen- 
tary sample of multijet events, where the azimuthal sepa- 
ration and -E™ ISS /m e ff, as well as all other event selection 
requirements, match those of the signal region, but where 
the r candidate is again randomly chosen from among the 
jet candidates. This provides an estimate of the multijet 
background yield in the signal region. It is found that the 
multijet background makes up only a few percent of the 
total SM background in the signal region. 

Possible contamination from SUSY signals has been 
considered in all three background-estimation control re- 
gions and is found to have a negligible effect on the results 
presented below. 

7. Systematic uncertainties 

Dominant systematic uncertainties on the estimated 
background yields are due to uncertainties in the jet en- 
ergy scale (3-8%) [7l|], jet energy resolution (6-13%) (zH, 
r energy scale (2-10%) [74|, statistical uncertainties in the 
data control regions (5-15%), and Monte Carlo uncertain- 
ties related to the extrapolation from the control regions 
to the signal region (10-20%). This last term includes sta- 
tistical uncertainties in the simulation, variations in the in 
the assumed W + jets/top/Z + jets mixture in the fake-r 
control region, and Monte Carlo generator uncertainties 
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(estimated by varying the shower matching, factorization 
and renormalization scales, a s , and the amount of initial- 
state and final-state radiation) [77j . Additional uncertain- 
ties on VF + jets and top-quark backgrounds are estimated 
by varying the assumed 6-quark identification efficiency 
within measured uncertainties (4-11%) [76|]. Uncertainties 
on the multijet background yield are estimated by study- 
ing correlations between m e g and the azimuthal separation 
between the leading two jets and the missing transverse 
momentum. Additional systematic uncertainties, includ- 
ing those on the pile-up description in the simulation, are 
considered and found to be negligible. 

In addition to the sources described above, systematic 
uncertainties on the SUSY signal cross section are esti- 
mated by varying the factorization and renormalization 
scales in PROSPINO up and down by a factor of two, by 
considering variations in as, and by varying the proton 
PDFs within their uncertainties. These theoretical uncer- 
tainties total typically 8-12% across the relevant region of 
parameter space. Uncertainties are calculated separately 
for individual SUSY production processes. 

8. Results 



Figure [2] shows the distributions of E™ ss , pT^, and m e s 
for data with all selection requirements applied except for 
that on m e ff, along with the corresponding estimated back- 
grounds. The numbers of expected SM background events 
and the observed number of events after the m e g require- 
ment are shown in Table [3] The data agree with the back- 
ground expectation. 

Based on these results, limits are placed on contribu- 
tions beyond the SM to the signal region. With 11 events 
observed and 13.2 ± 4.2 expected, an upper limit of 8.5 
on the number of events observed due to non-SM sources 
is derived at 95% confidence level (C.L.). This limit cor- 
responds to an upper limit on the visible cross section of 
4.0 fb, where the visible cross section is defined as the 
product of production cross section, branching fraction to 
at least one r lepton, acceptance, and efficiency using the 
event selection defined in Section [SJ For the two bench- 
mark points A = 30, tan/3 = 20 and A = 40, tan ,3 = 30 
the product of branching ratio to r-leptons, the accep- 
tance and the efficiency for this selection amounts to 1.47% 
and 1.69%, respectively. Figure [3] shows an interpreta- 
tion of the result as a 95% C.L. exclusion limit in the 
M mess = 250 TeV, N 5 = 3, /i > 0, C grav = 1 slice of the 
GMSB model. Figure [3] also shows the variation of the 
expected limit in response to ±ler fluctuations in the ex- 
pected SM background and the SUSY cross sections. The 
excluded regions are calculated using a profile likelihood 
method with systematic uncertainties modelled as varying 
Gaussian-distributed nuisance parameters [78l [79|. The 
resulting limit is compared with previous exclusion limits 
from searches for f and e production and GMSB topologies 
at LEP. The region of small A and large tan (3 is theoret- 
ically excluded since it leads to tachyonic states. In this 



model, the production of supersymmetric particles can be 
excluded at 95% C.L. up to A = 30 TeV, independent of 
tan/3, and up to A = 43 TeV for large values of tan/3. 
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Figure 2: Distributions of E™ lss , p^,, and m ff for data with all selec- 
tion requirements except for that on m e g , along with the correspond- 
ing estimated backgrounds. Backgrounds are taken from simulation 
and normalised with control regions in data. The solid (red) line 
with shaded (yellow) error band corresponds to the total SM predic- 
tion, while the points are data. The error bands indicate the size 
of the total (statistical and systematic) uncertainty. The notation 
GMSB(40,30) stands for the GMSB model with A = 40 TeV and 
tan^ = 30 and analogously for GMSB(30,20). 



9. Conclusions 

In conclusion, this Letter presents a search for super- 
symmetry in final states containing jets, missing transverse 
momentum, and at least one r lepton with the ATLAS 
experiment in ^/i = 7 TeV proton-proton collisions at the 
LHC. This is the first search in these final states at the 
LHC that includes events with one r lepton. No excess of 
events is seen beyond the expected Standard Model back- 
grounds in 2.05 fb -1 of data. Limits are placed on the 
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Table 3: Expected SM background event yields and number of events observed in data after the final requirement on m e ff. All systematic 
uncertainties are included here, and the uncertainty on EsMi the sum of all SM backgrounds, takes correlations between the individual 
background uncertainties into account. The truc-r purities are 53% and 64% for the top and W + jets backgrounds, respectively, and are 
negligible for the Z+jets and multijet backgrounds. For comparison, the estimated event yield for a GMSB signal with A = 40 TeV, tan /3 = 30 
is an additional 9.1 ± 1.7 events. 



Top TU+jets Z+jets Multijet £ SM Data 
5.6 ± 1.4 4.7 ± 1.5 2.4 ±0.7 0.5 ± 0.6 13.2 ±4.2 TT~ 
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Figure 3: Expected and observed 95% C.L. exclusion limits in the 
M mess = 250 TeV, N 5 = 3, n > 0, C grav = 1 slice of GMSB, together 
with the most stringent previous limits from OPAL [32| . The identity 
of the NLSP is indicated, with CoNLSP the region where the f and 
I are nearly degenerate. 

visible cross section and in the context of GMSB mod- 
els. The limits obtained extend the results from previous 
experiments. 
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